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A method for putting chiral fermions on the lattice 
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We describe a method to put chiral gauge theories on the lattice. Our method makes heavy use of the effective 
action for chiral fermions in the continuum, which is in general complex. As an example we discuss the chiral 
Schwinger model. 
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1. INTRODUCTION 

In spite of the well-known importance of chiral 
gauge theories there is very little nonperturbative 
information available about such models. Stan- 
dard nonperturbative techniques seem to fail, in 
particular the lattice approach faces great diffi- 
culties due to the notorious doubling problem. 
To overcome these difficulties it has been sug- 
gested to keep the fermionic degrees of freedom 
in the continuum and to latticize only the gauge 
fields ||l| . However, the details of such a procedure 
have not yet been worked out. In the present note 
we want to outline a few steps along this line (see 
also @). 

2. EFFECTIVE ACTION 

First of all we need the effective action for chi- 
ral fermions in the background of a continuum 
gauge field. This has been discussed by several 
groups of authors |^-|^. We shall mainly follow 
the approach of Alvarez-Gaume et al. 

Working in euclidean space we choose hermi- 
tian 7-matrices, 



and define 75 such that 

{75, 7m} = O.7I = 1 ,7-5^ = 75 



(1) 



(2) 



So the chiral projectors are given by 

^± = ^(1±75)- (3) 

The gauge fields are taken to be antihermitian: 

A;{x)^-A,{x). (4) 

We write the chiral Dirac operator as 

iDiA)^i{^+4P+) . (5) 

Consequently, both chiralities of fermions are 
present, but the gauge field couples only to one of 
them. For the vector Dirac operator we use the 
standard notation 



i^(^)=i(^+/) 



(6) 



Formally, i.e. up to regularization, the effective 
action W[A] can now be expressed in terms of the 
functional integral 



"D^jVip exp 



Aet[ib{A) 



d'^xi:{x)W{A)iJ;{x] 



(7) 



Since ^D{A) is not hermitian, W[A\ is expected to 
be complex in general. Up to local counterterms, 
the real part is given by the effective action for a 
Dirac fermion interacting with the gauge field A, 
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ReM/[A] = -ilndet(i^(A)) ^ 
and so it is gauge invariant. 



(8) 
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For the imaginary part one obtains 
IniM^[^] 

= l^'dTTr75 (^|-i^(A0) {ip{Ar))-' 

Here Tr denotes the operator trace, Ar inter- 
polates between Aq — and Ai = A (e.g. 
At — tA) and the current accounts for the dif- 
ference between the covariant and the consistent 
anomaly. Since we are ultimately only interested 
in anomaly free models, where contributions from 
these currents cancel, we shall omit the corre- 
sponding terms in the following. (Strictly speak- 
ing, the above formula can be derived only for 
perturbative gauge field configurations, i.e. con- 
figurations such that the topological charge van- 
ishes and the Dirac operator has no zero modes. 
However, one could more generally interpolate be- 
tween two configurations Aq and Ai to obtain 
ImM^[Ai] - ImW^[Ao].) Note that in eq. (|) only 
the vector Dirac operator appears. 

We can only mention in passing that ImT/F[j4] 
may be expressed in terms of a topological quan- 
tity, the spectral asymmetry or 77-invariant of a 
vectorlike Dirac operator in five dimensions (for 
a review, see [||). 

3. METHOD OF SIMULATION 

Since the effective action governing the dynam- 
ics of the gauge fields has turned out to be com- 
plex, one runs into well-known problems when at- 
tempting numerical simulations. At least in prin- 
ciple, however, one could incorporate ImVF[A] in 
the observable and update the gauge fields with 
the action 

Sg + RcW = 5*0 - I In dct(i 0) , (10) 

where So denotes the pure gauge field action. 
Remember that KeW is given by a vectorlike 
fermion action, whose lattice version should be 
unproblematic. 

Irrespective of the details of the simulation 
method one has to solve the question of how to 
calculate ImM^ for a given lattice gauge field con- 



figuration. A naive procedure would be to re- 
place p in eq. (||) by some lattice version of the 
Dirac operator, e.g. Wilson fermions. However, 
since ImW is related to a topological quantity, 
the r/-invariant, a look at the topological charge 
might be advisable. In that case such a proce- 
dure does not work: An expression for the lat- 
tice topological charge constructed out of a few 
plaquettes only suffers from large perturbative 
contributions leading to a divergent topological 
susceptibility in the continuum limit. The geo- 
metrical constructions of the lattice topological 
charge, on the other hand, which avoid this prob- 
lem, can be viewed as computing the continuum 
charge from a continuum gauge field derived from 
the lattice gauge field by a suitable interpolation 
Therefore we propose to define ImW by 
ImT^[y4] evaluated for this interpolated gauge po- 
tential A. A suitable interpolation has been ex- 
plicitly given in although the construction is 
not completely trivial as several constraints have 
to be fulfilled. In particular, gauge covariance has 
to be guaranteed. 

In order to really calculate ImVF[yl] one has to 
introduce some regularization. From the numeri- 
cal point of view it seems to be most convenient to 
use the lattice again, i.e. to compute on finer and 
finer sublattices of the original lattice anticipating 
that the results will converge. (The rigorous an- 
alytical investigations apply other regularization 
schemes, e.g. Pauli-Villars.) 

So one would proceed as follows. Start from 
a lattice gauge field U{x,iJ,) on the original lat- 
tice, whose spacing is put equal to one. Con- 
struct from U an interpolated continuum gauge 
field A. This is then used to calculate parallel 
transporters U {x, jj) on a finer lattice of spacing 
e < 1. These link matrices are raised to the power 
T to arrive at a lattice version of the gauge po- 
tential At ■ (An alternative consists in computing 
the parallel transporters on the finer lattice im- 
mediately from At-) As a regularized version of 
p{Ar)~^ we can then take the propagator of Wil- 
son fermions in the gauge field given by 

G{U^\x,y)^e-Hl'\U^\x,y) (11) 

with the fermion matrix 

M{U\x,y) 
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2e 
4r 



(12) 



Finally one has to calculate 



X lnU{x,fi)U{x,fiyG(U^\x + ft,x) 
1 

'2^" 



X Uix,fi)-^G{U^\x,x + fi) 



(13) 



in the limit e ^ 0, where tr denotes the trace over 
Dirac and internal indices. 

4. THE CHIRAL SCHWINGER MODEL 

As an example consider the chiral Schwinger 
model, i.e. a chiral U(l) model in two dimensions. 
On a square of side length L the fermionic part 
of the continuum action reads 



d^xij{x) (i ^~4{x)P+)tP{x). 



(14) 



where now 75 — 17172 and the potential Af^ {x) is 
real. We assume periodic boundary conditions for 
the gauge field, whereas the fermion field is taken 
antiperiodic. Choosing the linear interpolation 
Ar = tA so that 



ip{Ar) = i ^-t4{x) 



(15) 



we get for the relevant contribution to ImPF (up 
to regularization) 

|^'rfrTr75 

= '-j\TTvj,4{ip{Ar))-' 



1 



= |/ dT I d^xtTj5 4{x)G{Ar\x,x), 



'2™^-- ^ I- (16) 

^0 

where the propagator G{A\x, y) is determined by 
the equation 



I/)M)G{A\x,y)=8{x-y). 



(17) 



Since in two dimensions an explicit expression 
for G{A\x, y) is available we can evaluate (|l^ ) 
further. Being only interested in the gauge invari- 
ant contributions we fix the Landau gauge, 

dMx)^Q- (18) 

The zero momentum mode of the gauge field re- 
quires a special treatment, so we define 



fl/i ■= j2 J d^x^t^i^) 



(19) 



and assume —ixjL < < n/L. Introducing 
point splitting we find 

L 

d\ tr 75 4.{x)G{Ar \x,x + 5) 



X exp < — 1 



dz^A^{z)T 



= — {a2ni - ain2) / 

n ^ — ' 

X sin(L(6 + Ta) • n)e-"'-^'/'* 



47r \ - 

— (ain2 -I- jfli 



a2ni 



ia2) 



X e-(2Wi+&+^a)^/^27m/L + b + raf 

+ iL\aid2 - a2Si)/{7rS^) + Oi\S\) , (20) 

where 61 =62 ~ tt/L. Since the divergent term 
is odd in S, symmetric point splitting leads to a 
finite limit as (5 — > 0. Only the zero momentum 
mode of the gauge field contributes in Landau 
gauge. Note that working in infinite volume one 
does not have this mode and hence finds ImH^ = 
in Landau gauge (see, e.g., |^). 

Which results does our lattice recipe produce in 
this case? Unfortunately, an explicit expression 
for the lattice fermion propagator in a general 
background is not available. So we consider the 
special case of the gauge field configuration 



U{x,n) 



—ir/L < < TT / L . 



(21) 
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on the original lattice with lattice spacing 1. It 
should not come as a surprise that the inter- 
polated continuum gauge field turns out to be 



So we find for the parallel trans- 



porters on the sublattice of spacing e 

U{x,fi) =e'"''^ (22) 

and consequently 

In ij{x, n) = iea^ , U{x, ^i)^ = e'^'°f . (23) 

The required Wilson fcrmion propagator can eas- 
ily be calculated. According to our proposal we 
then have to compute 

~'^'^'^^{^^ - r) lnC7(x,/i) 

X U{x,fiyG{U^\x + fj.,x) 

+ tr75^(7^ +r) \nU{x,^i) 

xij{x,fi)-^G{if^\x,x + fi)} 

= 2e^^ [a2 cos(e/c2 + Tea2) sin(efci -I- rear) 

k 

— ai cos(efci -I- reai) sin(efc2 + ^£02)] 



E 

■ p 



(cos(efc^ -I- rea^) - 1) 



(24) 



where k runs over the momenta appropriate for 
our finite lattice. In the limit e — > one in- 
deed recovers the expression derived above by the 
point splitting technique in the continuum. So 
our recipe for the calculation of ImW from a lat- 
tice gauge field configuration works at least in this 
(almost trivial) case. 

5. OUTLOOK 

Obviously, there are many open problems 
which have to be solved before our proposal can 
become useful. For example, how can we deal 
with gauge field configurations leading to zero 
modes of the Dirac operator p{Ar) for some value 
of T? According to ref . ||^ , for every pair of eigen- 
values crossing zero, ImM^ picks up a contribu- 
tion ±7r. This problem has been ignored in our 



Schwinger model calculation, so there might be 
additional contributions to ImVF for certain con- 
figurations. If such configurations are statistically 
relevant, exp(i ImW) would fluctuate strongly 
and a Monte Carlo simulation might be hopeless, 
unless one invents a clever method to update with 
a complex action. In any case it would be inter- 
esting to see how other chiral fermion proposals 
deal with this problem. 

Even the calculation of ImW for a single config- 
uration poses several numerical challenges. E.g. 
one has to compute zero-mass fermion propaga- 
tors in a given background field. Furthermore, an 
operator trace has to be evaluated, which might 
be feasible with the help of a stochastic estimator. 

Let us close with the remark that the most el- 
egant approach would probably be to find a "ge- 
ometrical" expression for ImVF based on the re- 
lation to the 77-invariant. 
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